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ABSTRACT 

POINT- AGAPE is an Anglo-French collaboration which is employing the Isaac New- 
ton Telescope (INT) to conduct a pixel-lensing survey towards M31. Pixel lensing is a 
technique which permits the detection of microlcnsing against unresolved stellar fields. 
The survey aims to constrain the stellar population in M31 and the distribution and 
nature of massive compact halo objects (MACHOs) in both M31 and the Galaxy. 

In this paper we investigate what we can learn from pixel-lensing observables 
about the MACHO mass and fractional contribution in M31 and the Galaxy for the 
case of spherically-symmetric near-isothermal haloes. We employ detailed pixel-lensing 
simulations which include many of the factors which affect the observables, such as 
non-uniform sampling and signal-to-noise ratio degradation due to changing observ- 
ing conditions. For a maximum MACHO halo we predict an event rate in V of up 
to 100 per season for M31 and 40 per season for the Galaxy. However, the Einstein 
radius crossing time is generally not measurable and the observed full-width half- 
maximum duration provides only a weak tracer of lens mass. Nonetheless, we find 
that the near-far asymmetry in the spatial distribution of M31 MACHOs provides 
significant information on their mass and density contribution. We present a likeli- 
hood estimator for measuring the fractional contribution and mass of both M31 and 
Galaxy MACHOs which permits an unbiased determination to be made of MACHO 
parameters, even from data-sets strongly contaminated by variable stars. If M31 does 
not have a significant population of MACHOs in the mass range 10~ 3 Mq — 1 Mq 
strong limits will result from the first season of INT observations. Simulations based 
on currently favoured density and mass values indicate that, after three seasons, the 
M31 MACHO parameters should be constrained to within a factor four uncertainty in 
halo fraction and an order of magnitude uncertainty in mass (90% confidence). Inter- 
esting constraints on Galaxy MACHOs may also be possible. For a campaign lasting 
ten years, comparable to the lifetime of current LMC surveys, reliable estimates of 
MACHO parameters in both galaxies should be possible. 

Key words: dark matter — galaxies: haloes — galaxies: individual (M31) — Galaxy: 
halo — gravitational lensing. 



1 INTRODUCTION 



1.1 Conventional microlensing: landmarks and 
limitations 

The detection of the gravitational microlensing effect due to 
compact objects in the Galaxy is undoubtedly one of the 
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great success stories in astrophysics over the past decade. 
Surveys have discovered around 20 candidates towards the 
Magell anic clouds and several hundred towards the Galactic 
Bulge flUdalski et al. 1994 |Alcock et al. 1997j |Alard fc Guib- 



ert 199fJ Lasserre et al. 199£; Alcock et al. 2000). Amongst 
these candidates a number of exotic lensing phenomena have 
been catalogued, such as parallax effects, binary lensing 
(including spectacular examples of caustic-crossing events), 
and finite source-size effects. These discoveries are facilita ted 
by coo rdinated fo llow-up cam paigns such as PL ANET (Al- 
brow e | al. 1998| ) and MPS ( [Rhie et al. 1999| ) which act 



on microlensing alerts broadcast by the survey teams. The 
absence of certain microlensing signals has also yielded a 
clearer insight into the nature of halo dark matter. The null 
detection of short duration events towards the Large Mag- 
ellanic Cloud (LMC) by the EROS and MACHO surveys 
indicates that, for a range of plausible halo models, massive 
compact halo objects (MACHOs) within the mass interval 
10~ 7 — 10~ 3 M^ provide less than a quarter of the dark mat- 
ter (Alcock et al. 1998). This is an important result when 
set against the current insensitivity of other techniques to 
this mass range. 

Despite these successes a number of unsolved problems 
remain. The optical depth measured towards the Galactic 
Bulge is at least a factor two larger than can be accom- 
modated by theoretical models (e.g. Bissantz et al. 1997; 
Sevenster et al. 1999). Towards the LMC the rate of de- 
tected events is consistent with the discovery of a signif- 
icant fraction of the halo dark matter. However, the im- 
plied lens mass range (0.1 — 1 Mq) is not easily reconciled 
with e xisting con straints on baryonic dark matter candi- 
dates (Carr 1994), though the MACHOs need not neces- 



sarily be baryonic. Furthermore, the discovery of two pos- 
sible binary caustic-crossing events towards the LMC and 
the Small Magellanic Cloud (SMC) has thrown into ques- 
tion the very existence of MACHOs. Their caustic-crossing 
timescales, which provide an indicator of their line-of-sight 
position, seem to exclude either as being of halo origin, a 
statistically unlikely occurrenc e if the halo comprises a sig- 
nificant MACHO component (Kerins & Evans 1999). As a 
result, there is a growing body of opinion that all events 
observed so far towards the LMC and SMC may reside in 
the clouds themselves. However, this explanation is itself 
problematic because it requires that the clouds must either 
have a higher MACHO fraction than the Galaxy or comprise 
substantial but diffuse stellar components not in hydrody- 
namical equilibrium (Evans & Kerins 2000, and references 
therein). 

These problems highlight two principal constraints on 
the ability of conventional microlensing experiments to de- 
termine the nature and distribution of MACHOs in the halo. 
The first limitation is their inefficiency in differentiating be- 
tween lensing by MACHOs and self-lensing by the source 
population, since for most events one observes only a dura- 
tion and a position on the sky. These observables are only 
weakly correlated with the location of the events along the 
line of sight. The second constraint is the limited number 
of suitable lines of sight through the halo. Conventional mi- 
crolensing surveys require rich yet resolved stellar fields and 
are thus limited to just two lines of sight, the LMC and 
SMC, with which to probe MACHOs. The line of sight to 
the Galactic Bulge is dominated by bulge and disc lensing. 
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Figure 1. The principle of near-far asymmetry. The optical depth 
through the halo towards the far disc is larger than towards 
the near disc owing to the tilt of the disc confined within the 
spheroidal distribution of MACHOs. The distribution of Galaxy 
MACHOs, disc self-lensing events and variable stars does not ex- 
hibit asymmetry. 



The paucity of halo lines of sight, together with the rather 
weak dynamical and kinematical constraints on Galactic 
halo structure, also diminishes the prospect of being able 
to decouple information on the Galactic distribution func- 
tion and MACHO mass function. 



1.2 Beyond the Galaxy: a new target, a new 
technique 

The possibility of detecting MACHOs in an external galaxy, 
specifically M31, was initially explored by Crotts (1992) 
and by Baillon et al. (1993). Crotts (1992) pointed out 
that the high inclination of the disc of M31 would result 
in an asymmetry in the observed rate of microlensing if 
the disc is surrounded by a MACHO halo, as illustrated 
in Figure |l[ The fact that the M31 MACHO microlensing 
rate should be lower towards the near side of the disc than 
the far side, which lies behind a larger halo column den- 
sity, means that the presence of MACHOs in M31 can be 
established unambiguously. In particular, neither variable 
stars nor stellar self-lensing events in the disc of M31 should 
exhibit near-far asymmetry. Additionally, the external van- 
tage point serves to reduce systematic model uncertainties 
in two ways. Firstly, it permits a more accurate determina- 
tion of the rotation curve and surface brightness profile than 
is possible for the Galaxy, which reduces the prior param- 
eter space of viable galactic models. Secondly, it provides 
many independent lines of sight through the halo of M31, 
allowing the MACHO distribution across the face of the disc 
to be mapped and thus the halo distribution function to be 
constrained more or less directly. 

As pointed out by Baillon et al. (1993), another appeal 
of directing observations towards more distant large galaxies 
like M31 is the increase in the number of potential source 
stars, more than a factor of one thousand over the number 
available in the LMC and SMC, and all confined to within a 
few square degrees. However, this also presents a fundamen- 
tal problem in that the source stars are resolved only whilst 
they are lensed (and even then only if the magnification is 
sufficiently large). The presence of many stars per detector 
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pixel means it is often impossible to identify which is being 
lensed. Furthermore, the flux contribution of the unlensed 
stars dilutes the observed flux variation due to microlensing. 
Nonetheless, Baillon et al. (1993) determined from numerical 
simulations that the number of observable events, due to ei- 
ther the lensing of bright stars or high magnification events, 
is expected to be large. As a result of these studies, the An- 
dromeda Galaxy Amplified Pixel Experiment (AGAPE) and 
another group, Colu mbia- VATT, comm e nced observing pro- 
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grams towards M31 ( Ansari et al. 
19970 

One of the biggest technical difficulties facing surveys 
which look for variable sources against unresolved stellar 
fields is how to distinguish between flux variations due to 
changing observing conditions and intrinsic variations due 
to microlensing or stellar variability. For example, changes 
in seeing induce variations in the detected flux within a pixel. 
One must also deal with the consequences of positional mis- 
alignment between exposures, spatial and temporal varia- 
tions in the point spread function (PSF) and photometric 
variations due to atmospheric transparency and variable sky 
background. 

AGAPE has employed the Pixel Method to cope wit h 
the changing observing conditions (Ansari et al. 1997). 



AGAPE thoroughly tested this technique with a three- year 
campaign using the 2m Bernard Lyot telesc o pe at Pic du 



Mid i from 1994 to 1996 ( [Ansari et al. 1997j ; [Ansari et al 
199S |; |fje Du 2000[ ). Six fields covering about 100 arcmin^ 
centred on the bulge of M31 were monitored. Whilst the field 
of view was insufficient to conclude much about the nature 
of MACHOs, 19 candidate events were detected, though it 
is still premature to rule out many of them being intrinsi- 
cally variable sources, such as Miras or novae. One event, 
AGAPE Zl, appears to be a convincing lensing candidate 
as its flux increase and colour are inconsistent with that of 



a Mira or nova (Ansari et al. 1999). A longer baseline is 



needed to determine how many of the other candidates are 
due to microlensing. 

A major observing programme began on the 2.5m Isaac 
Newton Telescope (INT) in La Palma in the Autumn of 
1999, with a run of one hour per night for almost sixty nights 
over six months. The POINT-AGAPE collaboration is a 
joint venture between UK-based astronomers and AGAPE 
(where POINT is an acronym for "Pixel-lensing Observa- 
tions with INT"). We are exploiting the 0.3 deg 2 field of 
view of the INT Wide-field Camera (WFC) to map the dis- 
tribution of microlensing events across a large region of the 
M31 disc. Our initial observations of M31 with the INT em- 
ployed a V filter and the simulations reported here have been 
undertaken with parameters appropriate to V-band obser- 
vations. The strategy employed for the actual M31 mon- 
itoring campaign involves observations in three bands, g, 
r, and i [very similar to the bands employed by SLOAN 
(Fukugita et al. 1996)]. The multi-colour observations will 
improve our ability to discriminate against variable stars and 
the gri-Slter plus CCD combination offers a significant im- 
provement in sensitivity (the (/-band zero-point is some 0.4 
magnitudes fainter than that for V) . The simulation param- 
eters are thus somewhat conservative in this regard. The pro- 
gramme is being conducted in consort with the Microlensing 
Exploration of the Galaxy and A ndromeda (MEGA) survey 
(Crotts, Uglesich fe Gyuk 199S), the successor program to 



Columbia- VATT. Whilst POINT-AGAPE and MEGA are 
sharing the data, different techniques are being employed 
to search for micro lensing event s. Henceforth we use the 
term pixel lensing (Gould 1996) to describe microlensing 
against unresolved stellar fields, regardless of the detection 
technique. 

Whilst the technical viability of pixel lensing is now 
clearly established, a number of important theoretical issues 
are still outstanding. The principal concern is that the main 
observable in classical microlensing, the Einstein crossing 
time, is generally not accessible in pixel lensing. The Ein- 
stein crossing time is directly related to the lens mass, its 
transverse velocity and the observer-lens-source geometry. 
In pixel lensing the observed timescale depends upon addi- 
tional factors, such as the local surface brightness and the 
source luminosity and magnification, so the dependence on 
lens parameters is much weaker than for classical microlens- 
ing. 

The first detailed study of pixel lensing was undertaken 
by Gould (1996). He defined two regimes: a semi-classical 
regime in which the source star dominates the pixel flux and 
the observable timescale provides a fair tracer of the Ein- 
stein crossing time; and the "spike" regime where only high- 
magnification events are identified, and the timescales are 
only weakly correlated with the underlying Einstein cross- 
ing duration. Remarkably, Gould showed that, despite the 
loss of timescale information, in the spike regime one can 
still measure the microlensing optical depth. Using Gould's 
formalism, Han (1996) provided the first pixel event rate es- 
timates for the M31 line of sight. However, Gould's formal- 
ism assumes a fixed sampling rate and unchanging observing 
conditions. As such it is of limited applicability to a ground- 
based observing program. Gondolo (1999) has proposed an 
optical depth estimator based on the observed pixel event 
timescale. Whilst this estimator can be readily employed by 
a ground-based campaign, it is somewhat sensitive to the 
shape of the source luminosity function and is valid only 
to the extent that this can be taken to be the same for 
all source components. More recently, Baltz & Silk (1999) 
derived expressions for the pixel rate and timescale distri- 
bution in terms of the observable timescale, rather than the 
Einstein crossing time. Again, their study assumes constant 
sampling and observing conditions, as would be the case for 
space-borne programmes. 

Whilst these studies provide a solid foundation for pre- 
dictions of pixel-lensing quantities (i.e. timescales, rates and 
optical depth), none of them address to what extent one 
can constrain galactic and lens parameters, in particular the 
MACHO mass, from pixel lens observables. Gyuk & Crotts 
(2000) have shown that a reliable measure of the optical 
depth from pixel lensing can be used to probe the core ra- 
dius and flattening of the M31 MACHO halo. 

In this paper we quantitatively assess the degree to 
which the POINT- AGAPE campaign directed towards M31 
will constrain the fractional contribution and mass of the 
MACHOs. Since the answer inevitably depends upon the as- 
sumed galactic distribution function, we focus attention here 
on the simple case of spherically-symmetric near-isothermal 
halo models. The line of sight towards M31 is sensitive to 
two MACHO populations, our own and that in M31 itself, 
so we investigate the extent to which they can be distin- 
guished and probed independently. We also model the ex- 
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pected background due to variable stars and lenses residing 
in the disc and bulge of M31. 

The plan of the paper is as follows. In Section ^ we sum- 
marize the basic principles of pixel lensing, with emphasis on 
the differences between pixel lensing and classical microlens- 
ing. We describe our Monte-Carlo pixel-lensing simulations 
in Section including our event selection criteria and the 
incorporation of realistic sampling and observing conditions. 
In Section ^ we construct a reference model for the lens and 
source populations in the halo of the Galaxy and the halo, 
disc and bulge of M31, seeking consistency with the observed 
M31 rotation curve and surface brightness profiles. In Sec- 
tion P we present predictions for the POINT- AGAPE survey 
based on our simulations. In Section ^ we use the simula- 
tions to generate artificial data-sets and we investigate to 
what extent the MACHO mass and fractional contribution 
in the two galaxies can be recovered from the data. The 
results are summarized and discussed in Section Q 



2 PRINCIPLES OF PIXEL LENSING 

We review here some of the main aspects of pixel lensing and 
its differences with classical microlensing. A more detailed 
overview can be found in Gould (1996). 



2.1 Detecting pixel events 

Whilst in classical microlensing one monitors individual 
sources, in pixel lensing the sources are resolved only whilst 
they are lensed. We can therefore only monitor the flux in 
each detector element rather than the flux from individual 
sources. If a star is magnified sufficiently due to a lens pass- 
ing close to its line of sight, then the total flux in the detector 
element containing the source star (due to the lensed star, 
other nearby unlensed stars and the sky background) will 
rise significantly above the noise level and be recorded as an 
event. 

Before treating seeing variations the sequence of images 
must be geometrically and photometrically aligned with re- 
spect to some reference image, 1Z, as described in Ansari et 
al. (1997). The variations remaining after alignment are pri- 
marily due to changes in seeing and source flux, including 
microlensing events. To minimize the effects of seeing we de- 
fine our base detector element to be a superpixel: a square 
array of pixels. A superpixel is defined for each pixel, with 
that pixel lying at the centre, so that neighbouring super- 
pixels overlap with an offset of one pixel. The optimal size 
for the superpixel array is set by the ratio of the size of the 
seeing disc on images obtained in poor seeing to the indi- 
vidual pixel size. The INT Wide-field Camera (WFC) has 
a pixel scale corresponding to 0'.'33 on the sky, whilst poor 
seeing at La Palma is ~ 2". Adopting a very conservative 
value of 2'.'4 for the worst seeing leads to an optimized choice 
of 7 x 7 pixels for the superpixel array. A larger array would 
overly dilute source variations, whilst a smaller array would 
be overly sensitive to changing observing conditions. 

Whilst seeing variations are reduced by binning the pho- 
ton count into superpixels, this by itself is not enough to 
make them negligible. Residual variations are minimized by 
the Pixel Method, in which a simple, empirically-derived sta- 
tistical correction is applied to each image to match it to the 



characteristics of the reference image 1Z. The Pixel Method 
is discussed in Ansari et al. (1997) and described in detail by 
Le Du (2000). The method strikes a good balance between 
computational efficiency and optimal signal-to-noise ratio, 
with the resulting noise level approaching the photon noise 
limit. 

After alignment and seeing corrections the excess super- 
pixel photon count A7V p i x on an image i obtained at epoch 
ti due to an ongoing microlensing event is 



AJVp ix (ti) = iV b lL4 pix (ti) - 1] = f se eN s [A(ti) - 1]. 



(1) 



Here N s and Nbi are the source and baseline photon counts 
in the absence of lensing, A is the source magnification fac- 
tor due to lensing and / acc is the fraction of the seeing disc 
contained within the superpixel. The baseline photon count, 
iV u = Nggi(R.) + N sky (TZ), is the sum of the local M31 sur- 
face brightness (including N s ) and sky background contri- 
butions on the reference image. Whilst the quantities iVbi 
and f scc N s (A — 1) can be determined independently, N s and 
A cannot in general be inferred separately. It is therefore 
convenient to define A p i^ as the superpixel count variation 
factor, which acts as the observable analogue of A. 
The superpixel noise on image i is 



cr, = iaax[a<T(ti),aiNpi X (ti) 1 ' 2 ], 
where 

N pix (U) = AJVpi^ti) + N sky (t t ) + iV ga i 



(2) 



(3) 



refers to the superpixel photon count on image i prior to 
correction and, similarly, N sky and iV ga i are the uncorrected 
sky background and galaxy surface brightness contributions. 
The threshold noise level <jt is determined by the superpixel 
flux stability, and the scaling factor on takes account of the 
fact that the Pixel Method is not photon-noise limited. A 
preliminary analysis of a sequence of INT WFC images taken 
i n 1998 demon strated a flux stability level of 0.1 — 0.3% 
(Melchior 1999). We therefore adopt a conservative mini- 



mum noise level of <tt = 2.5 x 10 iVbi for our simulations. 
We also apply a constant scaling factor oti = 1.2, which is a 
l ittle larger t han typical for the AGAPE Pic du Midi data 
(Le Du 2000). In reality on varies slightly between images 
though we neglect this variation in our simulations. 

Note that iVg a i in equation ([]) is constant, despite the 
changing observing conditions. Though some variable frac- 
tion of the local patch of surface brightness is dispersed 
over neighbouring superpixels, the same amount of sur- 
face brightness leaks into the superpixel from neighbouring 
patches, so there is no net variation. The variation in N B k y 
results from changing moonlight and atmospheric trans- 
parency. 

We regard a signal as being statistically significant if 
it occurs at a level 3<7i above the baseline count N\>\. Our 
estimate of A?bi must be obtained from a sequence of im- 
ages and operationally is defined to be the minimum of a 
sliding average of superpixel photon counts over ten consec- 
utive epochs. In order for a signal to be detected on image 
i we therefore require a superpixel count variation factor 
^4pix(ti) > 1 + 3<Ti/7Vbi. From equation (|l]), a microlensed 
source satisfies this inequality provided that it is magnified 
by a factor exceeding 

;W*o = i + 7 ^ r . (4) 
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A special case of equation (|4j 
a threshold magnification of 

iVbi 



A T = 1 + 0.0075 



occurs when ai = err, giving 



(•») 



The extent to which residual temporal variations in / see and 
iVbi remain after image processing determines the factor by 
which ai exceeds the photon noise limit, so this excess noise 
is explicitly accounted for in equation m). 

Equation (Q) illustrates some important characteristics 
of pixel lensing. Firstly, pixel lensing does not depend di- 
rectly on the local surface brightness or sky background, 
but it does depend on their contribution to the noise at. 
Secondly, if the exposure time T exp is short, or the source 
star constitutes only a small fraction of the superpixel flux, 
so that N B -C cri, only rare high-magnification events are 
detected. The relationship between lens magnification and 
lens-source impact distance (measured in the lens plane) is 
as for the classical case: 



A = 



iVu' 2 + 4 



(6) 



where u is the impact distance in units of the Einstein radius. 
The maximum value for the impact distance can be obtained 
by inverting equation ([]) for A = Amia' 



,1/2 



An 



2 

min 



1/2 



(A min £ 10). (7) 



For pixel lensing in M31 we are often in the regime where 
N s <C cri because the source flux is much less than that of 
the galaxy and background, so it is not unusual to require 
A m in <; 10. In this case equations (^) and (Q) imply 



fsccN s / S eeiV; 
3<Ti 



3N, 



1/2 



(A min £ 10), 



(8) 



Since u max -C 1 [typically u max ~ 0(1O~ 2 — 10~ 3 )] only a 
small fraction of classical (u < 1) microlensing events are 
detectable. 

The dependence of it max on N s means that the pixel 
event rate depends on the source luminosity function 4>{M), 
the number density of sources in the absolute magnitude 
interval (M, M + dM). We can compute a theoretical upper 
limit, F p , for the pixel-lensing rate at sky coordinate (x,y) 
by taking A m i n = At so that it max = ii(At) = ut- In this 
case 



(9) 



where x and y axe Cartesian coordinates centred on M31 
and aligned respectively along the major and minor axes 
of the projected light profile. We define y to be positive 
towards the near side of the disc. The quantity F c is the 
classical (u < 1) event rate integrated over lens and source 
populations "( [driest 1991 ; Kiraga & Paczyfiski 1994), and 

Ju T (M,x,y)<t>(M)dM 
<UT(X ' y)) ^ j m)d M (10) 

is the mean threshold impact parameter at (x, y) averaged 
over </>. 

Whilst useful in providing a rough order of magnitude 
estimate, T p cannot be compared directly with observations 
because it assumes perfect sensitivity to all event durations 



and it also assumes that observing conditions are unchang- 
ing. Since one usually has j4 m i n > At, equation ( p"(i| ) also 
tends to overestimate the true mean pixel-lensing cross- 
section. One can regard r p , evaluated under the best ob- 
serving conditions, as providing a strict theoretical upper 
limit to the observed event rate, in much the same way as 
T c provides an upper limit to the observed rate in classical 
lensing. In Section ^ we set about obtaining a more realistic 
estimate of the observed pixel lensing rate. 



2.2 Degenerate and non-degenerate regimes 

In classical microlensing the most important observable is 
the Einstein radius crossing time, since this is directly re- 
lated to the position, motion and mass of the lens. Can we 
obtain similar information from the duration of pixel events? 

For a lens moving at constant velocity across the line of 
sight, u evolves with time t as in the classical case: 



u(tf = u(t ) 2 + ( 



t-to 



(11) 



where to is the epoch of minimum impact distance and t e is 
the Einstein radius crossing time. From equations (|^) and 
(p"l|), t c gives the timescale over which the source magnifica- 
tion A varies significantly. For large magnifications u ~ A~ 
from equation (Q), and inserting equation ([n]) into equa- 
tion (hi) gives 



AiV pix (i) 



fsec -^rnax 



\A{t) ;> io], 



(12) 



t c A 



where j4 max = A(to) is the maximum magnification. We in- 
fer that in pixel lensing the timescale over which the signal 
varies significantly is t e j4ma, x rather than t c . This means that, 
in the high-magnification regime, the pixel-lensing timescale 
bears little relation to t . We also see that the light-curve 
is degenerate under transformations A max — » Q^ m ax , 
N s /a and t c — > at c ( Wozniak & Paczynski 1997). So neither 
nor iV a can be determined independently. It may 
sometimes be possible to brea k this degeneracy b y looking at 
the wings of the light-curve (Baltz & Silk 199£), where dif- 



ferences between the true magnification and its degenerate 
form can become apparent. From equation ([]), the differ- 
ence between the exact expression for A{u) — 1 appearing in 
equation Qlf) and its degenerate approximation, is 



A(A - 



u + 2 

Uyju 2 + 4 



1 



■ill 
~~8 



(u<,l). (13) 



To discriminate reliably (say at the 3 a level) between the 
degenerate and non-degenerate cases requires / see iV s |A(A — 
1)| > 3<7j, so for the high- magnification regime we can write 
the condition for non-degeneracy as 



(«< 1). 



(14) 



Equation ([w]) demands that the superpixel noise be no 
greater than the contribution of the unlensed source to the 
superpixel flux. In general this will not be the case, so obser- 
vations will not be able to break the light-curve degeneracy 
and thus will not directly probe the Einstein crossing time. 
Since the underlying duration t c is not generally 
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measurable we use the observed full-width half-maximum 
(FWHM) event duration: 

1/2 



£fwhm = 2y/2t c 



a + 2 



a + 1 



Va 2 + 4a Va 2 + 2a 



(15) 



where a = ^4 max — 1. Since yl max for detected events is typ- 
ically larger in regions of higher surface brightness, and for 
fainter stars, tFWHM is correlated both with the disc surface 
brightness and the source luminosity function. This means 
that it is less strongly correlated than t e with the lens mass 
and velocity and the lens and source distances. 

The observed duration, £fwhm , does not afford us with 
as direct a probe of lens parameters as i c . We are therefore 
forced to rely on other observables, such as spatial distribu- 
tion, in order to probe the underlying MACHO properties. 
For M31 MAC HOs one can test for near- far asymmetry in 
the event rate (|Crotts 1992]). For Galaxy MACHOs there 



is no comparable signature. Looking from the centre of the 
Galaxy towards M31 the halo density distribution in the 
two galaxies is highly symmetric about the observer-source 
midpoint. Since the microlensing geometry is also symmetric 
about the midpoint the timescale distributions for Galaxy 
and M31 MACHOs are similar for the same mass function. 
Since our displacement from the Galactic centre is only 8 kpc 
(small compared to the scale of the haloes and the Galaxy- 
M31 separation) this geometrical symmetry is largely pre- 
served at our location. However, the Galaxy MACHO dis- 
tribution ought to be less concentrated than that of stellar 
lenses. One might hope to see this as an excess of events 
at faint isophotes which remains the same towards both the 
near and far disc. If MACHOs exist, the overall pixel-lens 
distribution will be superposition of several lens populations 
(Galaxy halo, M31 halo, disc and bulge) together with vari- 
able stars which, at least in the short term, appear indistin- 
guishable from microlensing. The task of disentangling each 
is therefore potentially tricky. 



3 SIMULATING PIXEL EVENTS 

A straightforward method for probing the lens populations 
is to construct simulations of the expected distribution of 
events for a particular telescope configuration, set of observ- 
ing conditions and selection criteria and then compare these 
predictions to observations. To this end we have constructed 
a detailed simulation of a realistic pixel- lensing experiment. 

Our simulation works by first computing a theoretical 
upper limit to the pixel rate for assumed M31 and Galaxy 
models. This estimate provides the basis for generating trial 
pixel microlensing events for which light-curves are con- 
structed and selection criteria applied. The precise details 
of our input galaxy models are discussed in Section ^ in 
this section we lay down the general framework for the sim- 
ulation. For each generated trial event, a pixel light-curve is 
constructed using a realistic distribution of observing epochs 
interrupted by poor weather and scheduling constraints. The 
effects of the sky background and seeing are explicitly taken 
into account in computing flux realizations and errors for 
each "observation". The observing sequence is then exam- 
ined to see whether the event passes the detection criteria 
— if it does, then the trial counts as a detected event. The 
simulation is terminated once 10 4 events are detected or 



10 trials generated, whichever is reached first. The fraction 
of trial events which are detected is used to compute the 
observed pixel rate. The statistical error on the rate deter- 
mination is typically about 3%. 

3.1 Generating trial events 

As the starting point for our simulation we use the theoret- 
ical pixel event rate as a function of position, T p (x,y), de- 
fined by equation (||). This quantity, evaluated for the best 
seeing conditions, always provides an upper limit to the de- 
tection rate at a given location and is therefore convenient 
to use to generate trial events. We compute F Pt j over a grid 
of locations (x, y) for each combination j of lens and source 
population. Near the centre of M31, j = 1 ... 8 since there 
are two source populations (M31 disc and bulge) and four 
lens populations (Galaxy halo, M31 halo, M31 disc and M31 
bulge). Beyond 8 kpc the M31 bulge is not in evidence, so 
j = 1 . . . 3. 

Given the grid of T Pt j(x, y), one can write the probabil- 
ity of observing an event at location (x, y) as 



P(x,y) oc AxAy^Sj{x,y)T p ^{x,y), 



(16) 



where Sj is the source surface density at (a;, y) for lens- 
source configuration j, and Ax and Ay are the local x and y 
grid spacings (required only for non-uniform grids). P(x,y) 
therefore reflects the total event rate in a box of area Ax Ay 
centred on (x,y). The box should be sufficiently small that 
Sj(x, y) and F p j(x,y) provide good estimates of the source 
density and theoretical rate anywhere within it. Having fixed 
the event location, F p j is then used to select the lens and 
source components from the probability distribution 



P(j) = 



(17) 



Once the event location and lens and source popula- 
tions have been decided, the next choice is the line-of-sight 
distances to the lens, D%, and source, D B : 

P(D S ) oc p B (D s )D 3 a /2 / P(D 1 )dLh 



(18) 



P(Di) oc p\(Di)^D\(D s - A), 

where p\ and p B are respectively the lens and source mass 
densities. These distributions reflect the dependency of the 
microlensing rate F p ,j on D s , integrated over all possible D\, 
and on D\, for a given D s . Next we require the lens mass m 
and relative transverse speed Vt. The lens mass realization 
is generated from the distribution 



P(m) oc m 1//2 i/)(m), 



(19) 



since, in the absence of finite source-size effects, F p oc R e tp oc 
m 1 / 2 ?/;, where ip is the lens mass function (i.e. the num- 
ber density of lenses per unit mass interval) and Re is the 
Einstein radius. The transverse speed Vt(Vi,Ks) is drawn 
from the assumed velocity distributions Pi(Vi) and P S (V S ) 
(see section ^), with V\ and V s the lens and source three- 
dimensional velocity vectors. Since the microlensing rate F p 
is proportional to VtP\P s rather than just P\P B , each of our 
realizations must be weighted by Vt in computing the final 
detection rate. Finally, we also need to generate the source 
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absolute magnitude M (defined for some photometric band) . 
The dependency of F p on M derives from the luminosity 
function <f> and the threshold impact parameter ut, so we 
have 

P(M)oau T {M,x,y)4>{M). (20) 



3.2 Generating light-curves 

At this point we have only simulated events according to 
the underlying distributions which govern F p ; we have yet 
to take into account the distribution of observing epochs, 
variations in observing conditions, or candidate selection cri- 
teria. 

The observing season runs from the beginning of August 
to the end of January, so we adopt the duration of an observ- 
ing season to be AT = 180 days. We assume 60 scheduled 
observing epochs per season — approximately the number 
of nights awarded for our 1999/2000 season. To construct 
a realistic sequence of observing epochs we assume that the 
WFC is mounted on the telescope and available for two- week 
periods every four weeks and that, on average, 25% of sched- 
uled observations are precluded by bad weather. Periods of 
poor weather are superposed on our initial observing se- 
quence to obtain a final sequence which typically comprises 
40-50 epochs per season. In practice we expect to obtain ob- 
servations on more epochs than this, but for the purposes of 
these simulations we assume 40-50 as a conservative lower 
limit. For example during the 1999/2000 season we have had 
observations on 56 nights. 

The epoch of maximum magnification to and the min- 
imum impact parameter w(to) are both chosen at random. 
u(to) is selected from the interval [0, wt], where the thresh- 
old impact parameter ut is computed from equations (|s]) 
and (W) taking A m \ n — At- This is all that is required to 
generate the underlying microlensing light-curve. 

To compute the pixel light-curve, we must also model 
the galaxy surface brightness and sky background. The sim- 
ulations presented here are performed in the V band and 
we use the radially-averaged surface brightness profile in Ta- 
ble VI of Walterbos & Kennicutt (1987) to estimate the con- 
tribution to the pixel flux of the galaxy background at the 
event location. The assumed sky background corresponding 
to a dark sky is listed in Table [ll along with other INT de- 
tector and site characteristics. The sky background varies 
over lunar phase and we adopt a contribution to the sky 
background from the full moon equivalent to 10 3 tenth mag- 
nitude stars per deg 2 (c.f. Krisciunas & Schaefer 1991). The 
contribution is modulated according to the lunar phase. The 
lunar contribution to the sky background also depends upon 
whether the moon is above the horizon and on its angular 
distance from M31. Our assumed value is taken to be an 
average over the positional dependence, so the true varia- 
tion in the sky background will be somewhat larger than 
we consider. We also simplify the computation of the seeing 
fraction / aee by adopting a Gaussian PSF with a FWHM 
equal to the seeing of the reference image. The position of 
the PSF maximum for the reference image is selected at 
random within the central pixel of the superpixel array. 

Using our computed values for / soc , the INT detector 
and site characteristics summarized in Table [l], and the mi- 
crolensing parameters generated for each event, we construct 



Table 1. Adopted characteristics of the INT observing site 
and Wide-field Camera (WFC). The sky background is given in 
mag arcsec -2 and the superpixel dimension is quoted in pixels. 
The zero-point is given in terms of the apparent magnitude of a 
source which results in a 1 photon sec -1 detection rate. All mag- 
nitudes are for the V band. Our survey is now observing in g, r 
and i filters. For comparison, the sky background and zero-point 
in g are 22.2 and 26.0, respectively. 



Characteristic 



INT WFC 



Best seeing 0'.'8 

Worst seeing 2'.' 4 

Reference image seeing 1" 

Sky background 21.9 

Scheduled epochs per season 60 

Field dimensions 32' X 32' 

Pixel field of view 0'.'33 

Superpixel dimensions 7x7 

Zero-point 25.6 

Exposure time per field 760 sees 



superpixel light-curves via equation (gj). The error at each 
epoch i is given by equation (El) . Poisson realizations for the 
superpixel flux at each epoch are generated from N p i x (ti) 
and <r;. 



3.3 Selection criteria and the observed rate 

The adoption of selection criteria inevitably reduces the 
number of detected events, but they are necessary to mini- 
mize the number of contaminating non-microlensing signals. 
As in all microlensing experiments the selection criteria must 
be based upon the quality of the data and the characteristics 
of non-microlensing variations. Ultimately the criteria must 
be derived from the data themselves, so they are inevitably 
experiment-specific and evolve as the experiment progresses. 
For our simulations we impose criteria based l oosely on the 



previous AGAPE pixel-lensing at Pic du Midi (Ansari et al 
1997|; [Le Du 2000[) 



The principal criterion for the selection of microlensing 
events in our simulation is that one and only one signifi- 
cant bump be identified on the light-curve. The bump must 
comprise at least three consecutive measurements lying at 
least 3cr above the baseline superpixel flux. Quantitatively, 
the significance of a bump is defined by its likelihood 



i— j+n,n>3 

ibump = Yl P(e > e 1 \e>, > 3), 



(21) 



where O; = [-/V p j x (t;) — iVbi]/(Xj and P(Q) is the probability 
of observing a deviation at least as large as O by chance. 
For a Gaussian error distribution, P = |erfc(0/-\/2). Equa- 
tion (^T|) indicates that we evaluate P(@i) only when Oi > 3. 
For our simulations we demand that a candidate have one 
bump with — InLbump > 100 and no other bump with 
— In Lbump > 20. We further demand that the epoch of max- 
imum magnification to lies within an observing season; we 
reject candidates which attain their maximum brightness be- 
tween seasons, even if they last long enough for the tails of 
the light-curve to be evident. This helps to ensure a reliable 
estimate of the peak flux, and in turn the FWHM timescale 

iFWHM- 
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Table 2. Parameters adopted for the density and velocity distributions for components of the Galaxy and M31. The bulge model 
is adopted from Kent (1989). 









Rotation 


Velocity 


Component 


Mass normalization 


Scale lengths 


speed 


dispersion 


M31 bulge 


M b = 4 x 10 10 Mq 




30 km s _1 


100 km s" 1 




M/L B = 9 








M31 disc 


p d (0) = 0.2 M pc- 3 


H = 0.3 kpc, 


235 km s" 1 


30 km s" 1 




M/L B = 4 


h = 6.4 kpc 






M31 halo 


p h (0) = 0.23 M pc" 3 


a = 2 kpc, 





166 km s _1 






Rnmx = 200 kpc 






Galaxy halo 


p h (0) = 0.036 M pc" 3 


a = 5 kpc, 





156 km s _1 






iW = 100 kpc 







The bump criterion is both a signal-to-noise ratio con- 
dition and a test for non-periodicity. It is crucial for distin- 
guishing microlensing events from periodic variables, though 
long-period variables, such as Miras, may pass this test in 
the short term. In addition to the bump test, one can also 
test the goodness of fit of the light-curve to microlensing, 
which helps to distinguish microlensing from typical novae 
light-curves. Though the presence of the background means 
that pixel events will not in general be achromatic, the ra- 
tio of the flux increase to baseline flux in different colours 
should nonetheless be independent of time, so this provides 
another test for microlensing. Colour information may also 
help to exclude some long-period variables in the absence 
of a sufficient baseline of observations. In Section ^ we also 
exploit differences in spatial distribution to separate statis- 
tically lensing events from variable stars. 

For real data-sets we would require more criteria in or- 
der to avoid excessive contamination from variable stars. 
For now we are simulating only microlensing events, so we 
are assured of no contamination in our selection. However, 
the cuts adopted above would be responsible for many of 
the rejected candidates in a real experiment, so the absence 
of further criteria should not lead to a gross overestimate of 
the rate. In any case, we have been deliberately conservative 
with our choices of sky background level, worst seeing scale, 
the number of epochs per season and the pixel stability level 
ot- We therefore feel our predictions are more likely to be 
underestimates of the actual detection rate. 

The observed rate can be now readily computed from 
r p , the number of generated trials and the fraction of these 
whi ch p ass the detection criteria. As mentioned in Sec- 
tion 3.1, the way in which velocities are generated in the sim- 
ulations means that the correct rate is obtained by weight- 
ing each event by its transverse speed Vt- Thus, the observed 
rate for lens component j is 



^dot(j) 



v. 



CO 



v t , k 



(22) 



where {Tp,j)m,y is the spatial average of T Pi j (summed over 
source populations), the lower summation is over all Atrial 
trial events generated for lens component j and the upper 
summation is over the Ndet detected events which pass the 
selection criteria. The total number of events after n observ- 
ing seasons is 



N = n AT 10' 



0.4«M)-M gal ) V^pobs 
/ / PJ ' 

i 



(23) 



where (M) is the average absolute magnitude of the sources 
(integrated over the luminosity function) and M ga i is the 
absolute magnitude of M31 (M v = -21.2). 



3.4 Simulated light-curves 

Three light-curves generated for a first-season simulation in- 
volving 0.1 Mq MACHOs are shown in Figure [| The galac- 
tic models required for the simulation are discussed in Sec- 
tion ^j. The light-curves illustrate the range in signal-to-noise 
ratio. The down-time for the WFC is evidenced by the way 
in which the epochs are clumped into two- week periods. The 
variation in the size of the error bars reflects the simulated 
variation in observing conditions. 

Figure ^a shows an M31 halo lens magnifying a bulge 
star (My = —0.4) and is a typical example. The underly- 
ing maximum magnification for this event is A max = 18, 
whilst the maximum enhancement in superpixel flux is 
-Apix(to) — 1-06, indicating that the unlensed source is con- 
tributing less than 0.4% of the superpixel flux. For this event 
(fwhm = 5 days and t e = 28 days. Figure ^b, which il- 
lustrates a poor candidate with a low signal-to-noise ratio, 
involves a Galaxy MACHO and My = 1.8 bulge source 
contributing only 0.1% of the superpixel flux (yl max = 42, 
A p ix(to) — 1.05). In this example £fwhm = 5 days and 
t c = 68 days. Though there appears to be evidence of a 
second bump after the main peak these points are all within 
3<t of the baseline and so do not count as a bump. Fig- 
ure ^c shows a high signal-to-noise ratio "gold-plated" event 
in which a very luminous (My = —4) disc source is lensed 
by an M31 MACHO (A max = 5, A pix (t ) = 2.1) with an ob- 
served duration £fwhm = 19 days and underlying timescale 
t c = 33 days. Here the bright unlensed source accounts for 
27% of the superpixel flux. 



4 LENS AND SOURCE MODELS 

In order to make quantitative estimates for pixel-lensing ob- 
servables, we must specify models for the principal Galaxy 
and M31 lens and source components. For M31 the main 
populations are the bulge, the disc and the dark MACHO 
halo. For the Galaxy only the MACHO halo is important 
since the disc does not contribute significantly. Our complete 
model therefore consists of these four populations. Two pop- 
ulations, the M31 disc and bulge, also provide the sources, so 
in total we have eight different lens-source configurations. 
For each population we must specify distributions for the 
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(a) 



CO 
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a (arcmins) 
(b) 



o 




a (arcmins) 

Figure 3. (a) The overall surface brightness profile (solid line) as a function of semi-major axis a for our M31 model produced by 
the combined bulge (dashed line) and disc (dot-dashed line) light. The crosses are radially-averaged measurements from Table VI 
of Walterbos & Kennicutt (1987). (b) The overall rotation curve (solid line) for the same M31 model summed over bulge (dashed 
line), disc (dot-dashed line) and halo (dotted line) contributions. The crosses are from Figure 2 of Kent (1989) and are based on 
emission line measurements. For conversion to distance 1 kpc = 4.5 arcmin. 



density and velocity. Additionally, we must specify the lens 
mass and a luminosity function for the source populations. 
Throughout we assume a disc inclination of 77° and a dis- 
tance to M31 of 770 kpc, consistent with recent determina- 
tions (e.g. Stanek & Garnavich 1998). 

Whilst the present paper is concerned only with quan- 
tities relating to M31 and Galaxy MACHOs, we must 
nonetheless include other significant lens components in our 
modeling in order to properly characterize the complexity 
of extracting physical information from observations. For 
the observations, unlike the simulations, we do not know in 
which population a particular lens resides. 

The haloes are modeled as simple near-isothermal 
spheres with cores, having density profiles 



Ph 



Ph(0) 




— 7~ ', 1 

a + r 



(r < R n 
(r > Rr, 



(24) 



where ph(0) is the central density, a is the core radius, 7? ma x 
is the cutoff radius and r is the radial distance measured 
from the centre of either M31 or the Galaxy. The assumed 
values for ph(0), a and -R ma x are given in Table |2| The halo 
fraction determinations in Section [j] are made with respect 
to these density normalizations. In our model the M31 halo 
has about twice the mass of the Galactic halo, though this 
mass ratio is controversial and has been challenged recently 
by Evans & Wilkinson (2000) who have studied the kine- 
matics of several satellite galaxies around M31. 

The M31 disc is modeled by the sech-square law: 



pd = Pd(0) exp sech 2 (Jj^j , 



(25) 



where a is the radial distance measured in the disc plane and 
z is the height above the plane. The normalization Pd(0), 
scale-height H and scale-length h are given in Table Q. 

The bulge distribution is based on the work of Kent 
(1989). Kent models the bulge as a set of concentric oblate- 
spheroidal shells with axis ratios which vary as a function 
of semi-major axis. We use the tabulated spatial luminosity 
density values in Table 1 of Kent (1989) and normalize the 
bulge mass under the assumption that the light traces the 
mass (constant bulge mass-to-light ratio) . The mass normal- 
ization Mb is listed in Table |^. The assumption of axisymme- 
try may be over-simplistic since the misalignment between 
the disc and bulge position angles probably implies a triax- 
ial structure for the bulge. However, we are only indirectly 
concerned with bulge lensing in so much as it contaminates 
halo lensing statistics, so deviations from axisymmetry are 
not crucial. 

The rotation curve and surface brightness profile for the 
adopted M31 components are shown in Figure |j| In con- 
structing the surface brightness profile, we have assumed 
B-band mass-to-light ratios M/Lb = 4 for the disc and 
M/ Lb = 9 for the bulge, consistent with that expected 
for typical disc and bulge populations. The overall surface 
brightness profile is shown by the solid line in Figure |^a, with 
the disc and bulge contributions indicated by the dashed and 
dot-dashed lines, respectively. The crosses are the radially 
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averaged measurements from Table VI of Walterbos & Ken- 
nicutt (1987). In Figure the solid, dashed and dot-dashed 
lines show the overall, disc and bulge contributions to the 
rotation curve, with the dotted line giving the halo contri- 
bution. The crosses are from Figure 2 of Kent (1989) and are 
based on the emission- line curves of Brinks & Shane (1984) 
and Roberts, Whitehurst & Cram (1978). The fit to both 
the surface brightness and rotation profiles is good, given 
the simplicity of the models. 

The lens and source velocities are described by rota- 
tional and random components. The rotation velocity for 
each component is given in the 4th column of Table [| The 
random motions are modeled by an isotropic Gaussian dis- 
tribution with a one-dimensional velocity dispersion given 
by the 5th column. When calculating the relative transverse 
lens speed Vt> we take account of both the motion of the 
source and the observer. The observer is assumed to move in 
a circular orbit about the centre of the Galaxy with a speed 
of 220 km s" 1 . We do not assume any relative transverse 
bulk motion between the Galaxy and M31. In practice, only 
the observer's motion is of consequence for Galaxy lenses, 
and only the source motion for M31 lenses. 

Since one of the questions we wish to address is how 
well pixel-lensing observables can characterize the MACHO 
mass, we shall simply model the Galaxy and M31 MACHO 
mass distributions by a Dirac 5-function: 



ipOrrih) oc Sim — mh), 

m h 



(26) 



The stellar lens mass distribution in the disc and bulge is 
described by a broken power law: 



ip(m s 



m s " u -' b (mi < m s < 0.5 Mq) 

-™.-2.2 



(0.5 M < m s < m u ) 



(27) 



The mass function is normalized to yield the same value for 
^(0.5 Mq) for either slope. We take a lower mass cut-off 
mi = 0.08 Mq and an upper cut-off m u = 10 Mq, cor- 
respondin g closely to the local Solar ne ighbourhood mass 
function (|Gould, Bahcall fc Flynn 1997|). Whilst this is a 



reasonable assumption for stars in the M31 disc, the mass 
function will overestimate the contribution of massive stars 
in the older bulge. The higher M/ Lb assumed for the bulge 
also requires that the disc and bulge mass functions be dif- 
ferent. However, the slope at high masses is steep, so the 
contribution of high mass stars to the lensing rate is in any 
case small. Furthermore, as already mentioned, we are only 
interested in the bulge population as a contaminant of the 
halo lensing statistics. The choice of upper mass cut-off for 
the bulge is therefore not critical for the present study, so 
we simply adopt the same mass function for the disc and 
bulge. 

The stellar components provide both lenses and sources. 
We assume that the lens and source populations are the 
same and so described by the same density, velocity and 
mass distributions. For the disc and bulge sources, we use the 
1^-band luminosity function of Wielen, Jahreiss & Kriiger 
(1983) for stars with Mv > 5 and that of Bahcall & Soneira 
(1980) for Mv < 5. The two functions are normalized to the 
same value at Mv = 5. A more detailed study of the M31 
luminosity function is underway (Lastennet et al. 200C). 




MACHO mass (M ) 

Figure 4. The expected event rate as a function of MACHO mass 
for full MACHO haloes. The rates are averages over the M31 disc 
for M31 (solid line) and Galaxy (dashed line) MACHOs and are 
computed from ten seasons of data comprising 460 epochs. 



5 PREDICTIONS AND TRENDS FOR PIXEL 
LENSING 

The simulations for the POINT-AGAPE survey are per- 
formed over 1, 3 and 10 observing seasons for 9 MACHO 
masses spanning the range 10 -3 — 10 Mq. Each simulation 
produces an estimate of the number of events across the 
whole M31 disc for each lens component, together with a 
library of typically 10 4 candidates containing information 
such as the lens position, duration and transverse velocity. 
Since t e cannot generally be measured from the light-curve, 
we output both t c and £fwhm • The event libraries can be 
filtered to provide an estimate of the pixel-lensing rate for 
any field placement. 



5.1 Number of events 

Whilst the factor 10 3 gain over LMC/SMC searches in the 
number of sources certainly boosts the rate of events, the 
fact that M31 pixel-lensing searches can typically detect 
only high-magnification events means that the gain in the 
rate is not of the same order. Nonetheless, as Figure ^ indi- 
cates, the expected pixel-lensing rate is almost an order of 
magnitude larger than for current LMC/SMC experiments 
for same lens mass and halo fraction. In the figure we have 
plotted the expected number of events for M31 MACHOs 
(solid line) and Galaxy MACHOs (dashed line) per season 
per deg 2 , assuming MACHOs comprise all the halo dark 
matter of both galaxies. The rates are averages over the 
whole M31 disc (rather than for a specific field placement) 
determined from simulations spanning ten seasons and 460 
observing epochs. Within the first season the sensitivity to 
very massive MACHOs will be a little less than indicated in 
Figure ^J. 

The rate of events occurring within the two INT WFC 
fields for their first season (1999/2000) positions are dis- 
played in Table [| This excludes events occurring within 
5 arcmin of the centre of M31 because this region is domi- 
nated by stellar self-lensing (see Section 5.S). Only a couple 
of self-lensing events per season are expected outside the ex- 
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Table 3. The expected number of M31 and Galaxy MACHO 
detections per season (averaged over ten seasons comprising 460 
epochs) for a range of MACHO masses based on the placement of 
the two INT WFC fields in the 1999/2000 observing season. The 
numbers assume the haloes of both galaxies completely comprise 
MACHOs, though we exclude events occurring within 5 arcmin of 
the centre of M31. For comparison, the expected number of bulge 
and disc self-lensing events occurring outside the exclusion zone 
is 2.2 per season. The Monte-Carlo error for a given sequence of 
observing epochs is about 3%. 

Mass/M AT(M31)/yr N (Galaxy)/yr 
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Figure 5. Observed MACHO timescale distributions for a range 
of MACHO masses. The curves represent the combined M31 
and Galaxy MACHO normalized timescale distributions, shown 
in terms of the measured FWHM timescale tpwHM- From the 
lightest to the darkest curve the MACHO mass is 0.001 Mq, 
0.003 M Q , 0.01 M Q , 0.03 Mq, 0.1 Mq, 0.3 Mq, 1 M Q , 3 Mq 
and 10 Mq. 

elusion zone. The Monte-Carlo error in the values in Table ^ 
is small, only about 3%, but one should expect a larger vari- 
ation when comparing different seasons with different num- 
bers of epochs (in addition to Poisson variations). 

From Figure ^| and Table ^ we see that the sensitivity 
to MACHOs peaks at a mass around 0.003 — 0.01 Mq, when 
around 140 MACHO events can be expected within the INT 
WFC fields for full haloes. Below 10" 3 M finite -source size 
effects become important, so the expected number of events 
will drop off rapidly. At the high mass end, even haloes com- 
prising MACHOs as massive as 10 Mq provide a rate of 
several events per season. The number of M31 MACHOs is 
about twice as large as the number of Galaxy MACHOs for 
the same mass and fractional contribution, which is a direct 
consequence of the mass ratio of the halo models we adopt. 



Table 4. Average timescales for M31 and Galaxy MACHO pop- 
ulations for a range of masses. (tFWHlvl) is the mean FWHM 
duration, as measured from the light-curve, whereas (t e )det & n d 
(ie)pop are the mean Einstein radius crossing durations of de- 
tected events and of the underlying population, respectively. 



Mass/MQ 




M31 


Galaxy 


0.001 


(*fwhm) 


3.8 


4.0 




(te)det 


6.2 


6.3 




/+ \ 

\tc/pop 


Z.O 


O.l 


0.003 


(*fwhm) 


5.1 


5.1 




(*e)det 


9.1 


9.2 




It \ 


a n 


. O 


0.01 


(*fwhm) 


7.2 


7.8 




(*e)det 


14 


15 




It \ 


7 ^ 


7 


0.03 


(*fwhm) 


9.7 


9.4 




(ie)det 


21 


22 




It \ 




"1 7 


0.1 


(*fwhm) 


13 


13 




(*e)det 


34 


37 




It \ 


9^ 


OX 


0.3 


(*fwhm) 


16 


17 




(ie)det 


52 


57 




(ie)pop 


40 


53 


1 


(*fwhm) 


21 


23 




(*e)det 


82 


98 




(ie)pop 


73 


97 


3 


(^fwhm) 


26 


28 




(*e)det 


130 


160 




(ie)pop 


130 


170 


10 


(*fwhm) 


41 


32 




(*e)det 


220 


300 




(ie)pop 


230 


310 



5.2 Timescale distributions 

In Figure ^| we plot the timescale distributions for the de- 
tected MACHOs for a range of masses in terms of £fwhm- 
The distributions for nine MACHO masses, spanning four 
orders of magnitude, are plotted. The masses are as listed 
in Table ^, with darker lines corresponding to more massive 
MACHOs. Since the timescale distributions for Galaxy and 
M31 MACHOs are practically indistinguishable for a given 
mass, in Figure || we have combined their timescale distri- 
butions, so the normalization of each curve is determined by 
the combined pixel-lensing rate shown in Figure ^| for each 
halo. 

Whilst there is a clear trend of increasing *fwhm with 
increasing MACHO mass, the correlation is much weaker 
than for t c . For example, a duration £fwhm = 10 — 20 days 
is typical of a 0.1 Mq lens, but it is also not unusual for a 
lens as light as 10~ 3 Mq or as heavy as 10 Mq. Figure 
shows how the average duration (iFWHivi} varies with mass 
separately for M31 (solid line) and Galaxy (dashed line) 
MACHOs. Over four orders of magnitude in mass (£fwhm) 



© 0000 RAS, MNRAS 000, 000-000 



Theory of pixel lensing towards M31 I 13 




MACHO mass (M e ) 



Figure 6. The mean FWHM duration, (tpwHM)i as a function 
of MACHO mass. Line coding is as for Figure t&. 




MACHO mass (M Q ) 



Figure 7. The mean ratio of FWHM duration to Einstein radius 
crossing duration, (fpwHM/i:)i as a function of MACHO mass. 
Line coding is as for Figure W 

varies by about one order of magnitude, increasing from 
4 days for 10" 3 M MACHOs to 35 days for 10 M MA- 
CHOs (see also Table ^). For our sampling strategy we find 
empirically that (£fwhm) oc m^ 4 , whereas the average Ein- 
stein radius crossing timescale for the underlying population 
of microlensing events (with u < 1) scales as (t c ) P op oc "rnV 2 . 

The mean ratio (iFWHM/i e ) is displayed in Figure m for 
detected events. It is clear that the ratio is not fixed but 
steadily decreases with MACHO mass. For low MACHO 
masses with short durations, sampling imposes a lower limit 
on tFWHM and a loose lower limit on t c as well. Whilst most 
events involving ~ 1CF 3 Mq lenses are too short to be de- 
tected, those that are either have an unusually long t e or 
occur in regions of low surface brightness (which maximizes 
£fwhm for a given magnification). Thus (£fwhm /^e) is typ- 
ically larger for the observed events. At the other end of 
the mass scale the converse is true. The total observation 
baseline imposes a maximum cutoff in Jfwhm and a loose 
upper limit in t c . Those events which are detected either 
have an unusually short t or else tend to occur in regions 



of high surface brightness where £fwhm is minimized for a 
given magnification. So (tFWHM/ic) tends to be smaller for 
observed events. From Table ^ we see that the average dura- 
tion of detected events (t e )det does not trace the population 
average {t c ) pop . This is a consequence of sampling bias. 

5.3 Spatial distributions 

Since event timescales give only limited information in pixel 
lensing, the location of each event on the sky is a crucial 
observable. A robust measurement of near-far asymmetry 
in the event distribution would indicate the existence of an 
extended spheroidal population of lenses within which the 
visible M31 disc and bulge are embedded. Thus it would 
represent very firm evidence for the existence of MACHOs. 

In Figure |^ we display the distribution of events across 
the face of the M31 disc after three observing seasons for the 
case where the haloes of both M31 and the Galaxy are full 
of 0.3 Mq MACHOs. The axes are labeled in arcmins and 
are aligned along the major and minor axes of the disc light 
profile. The dashed-line templates indicate the positions of 
the two INT WFC fields for the 1999/2000 observing season. 

In Figure [s^, the positions of all detectable events are 
shown. MACHOs from the Galaxy halo are shown in green 
whilst M31 MACHOs are shown in blue. We find that within 
the central 5 arcmins (denoted by the circle) most events are 
produced by ordinary stellar lenses in the disc and bulge 
(shown in red). In Section |(| where we try to estimate 
MACHO parameters from simulated data-sets, we disregard 
events occurring within this region so as to minimize con- 
tamination from stellar lenses. 

Figure [|j shows only the M31 MACHO distribution. 
The excess of events between y = — 10 and —20 arcmins 
(along the minor axis towards the far side of the disc) com- 
pared to the number between y — +10 and +20 arcmins is a 
consequence of near-far asymmetry in the pixel-lensing rate. 
The strength of this asymmetry depends upon the number 
of M31 MACHOs which, in turn, depends upon their mass 
and density contribution, as well as the span of the observa- 
tion baseline. The presence of Galaxy MACHOs makes the 
asymmetry harder to detect, so the ratio of M31 to Galaxy 
MACHOs is another factor which determines whether or not 
the asymmetry is measurable. It is evident from the figure 
that there are very few events at \y\ J; 25 arcmin. This is 
due to the decrease in both the number of sources and the 
signal-to-noise ratio (because the sky background provides 
a larger fraction of the total superpixel flux). The presence 
of the sky background effectively imposes a cut-off in the 
spatial distribution. 

Figure ^| shows the spatial distribution for a range of 
MACHO masses expected after three seasons. We again as- 
sume that the MACHO mass is the same in both galaxies 
and that MACHOs provide all the dark matter in the two 
haloes. FigureHa is for a MACHO mass of 0.1 Mq. In Fig- 
ures §d and |c the MACHO mass is 1 Mq and 10 Mq 
respectively. The most obvious trend in the MACHO distri- 
butions is the decrease in the number of detectable events 
for models with more massive MACHOs. However, even for 
a mass as large as 10 Mq we still expect to detect 30 — 40 
MACHOs within the INT fields if they make up all the dark 
matter. After three seasons even these massive MACHOs 
out-number the disc and bulge lenses lying outside of our 
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Figure 8. A realization for the spatial distribution of pixel-lensing events after three seasons of observing, assuming MACHOs 
have a mass of 0.3 Mq and provide all the halo dark matter in the Galaxy and M31. The axis labelling is in arcmins. (a) The 
distribution of all events. The green dots represent the foreground Galaxy MACHO distribution, the red dots represent stellar lens 
events and the blue dots depict the M31 MACHO distribution. The circle centred on the origin demarcates the exclusion zone for 
the MACHO analysis, inside which the rate is dominated by stellar self-lensing. The dashed-line templates show the positions of 
the two INT fields for the 1999/2000 observing season, (b) The distribution of M31 MACHOs only. The near-far asymmetry can 
be seen by comparing event number densities at ±(10 — 20) arcmins along the minor axis. 



exclusion zone. This highlights one of the benefits of pixel 
lensing: the reduction in (fwhm due to the presence of many 
neighbouring unresolved sources means that more events 
with relatively large i e can be detected and characterized 
within a given observing period. In this respect, pixel lensing 
is relatively more sensitive to massive MACHOs than con- 
ventional microlensing experiments, which require resolved 
sources. 

Another noticeable trend in Figure ^| is that more mas- 
sive MACHOs are concentrated towards the central regions 
of the M31 disc. The main reason is that the MACHO and 
source surface densities are largest in this region, so the 
probability of an event occurring there is larger. However, 
another factor is that it is in the regions of highest surface 
brightness that the ratio £fwhm /te is minimized for a given 
magnification. For the 10 Mq MACHO model, where many 
events may have a duration t e exceeding the survey lifetime, 
this means more light-curves can be fully characterized, en- 
abling these events to be flagged as microlensing candidates 
within the observing period. The converse is true for low- 
mass MACHOs with short t e . Their distribution is biased 
towards regions of lower surface brightness where tFWHM /te 
is maximized. This effect provides a further degree of dis- 
crimination for different lens masses and means that, for 
example, a halo with a modest contribution of low mass 
MACHOs may be distinguished from one with a substantial 



fraction of more massive lenses, even if the number of events 
for the two models is comparable. This in part makes up for 
the fact that £fwhm is a less powerful discriminant than t B . 



6 ESTIMATING MACHO PARAMETERS 

In the previous section we found that, whilst the timescale 
information in pixel-lensing studies is somewhat more re- 
stricted than in conventional microlensing we do, at least for 
M31, have important information from the spatial distribu- 
tion of lenses. We now address to what extent pixel-lensing 
observables permit a reconstruction of the MACHO mass 
and halo fraction in the Galaxy and M31. 

6.1 Maximum-likelihood estimation 

Alcock et al. (1996) presented a Bayesian maximum likeli- 
hood technique to estimate the Galaxy MACHO mass and 
halo fraction from the observed event timescales towards the 
LMC. Evans & Kerins (2000) extended this to exploit the 
spatial distribution of observed events, and also to allow for 
more than one significant lens population. For pixel lensing 
towards M31 we must also consider the effect of contam- 
ination by variable stars. This is likely to be a significant 
problem in the short term. A baseline of more than three 
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Figure 9. Realizations for the spatial distribution of pixel-lensing events after three seasons of observing, assuming the MACHOs 
in M31 and the Galaxy have the same mass and provide all the halo dark matter, (a) The distribution for 0.1 Mq MACHOs; (b) 
1 M Q MACHOs; (c) 10 M Q MACHOs. The lines and symbols are as for Figure |. 



years should be sufficient to exclude periodic variables, such 
as Miras, but there still remains the possibility that, oc- 
casionally, the signal-to-noise ratio may be insufficient to 
distinguish between novae and microlensing events. By tak- 
ing account of variable stars in our likelihood estimator we 
allow ourselves to make an estimate of the MACHO mass 
and lens fraction which, even in the short term, is robust 
and unbiased. 

In order to allow for different MACHO parameters in 
the two galaxies we propose an estimator which is sensitive 
to five parameters: the MACHO mass and halo fraction in 
both the Galaxy and M31, and the degree of contamination 
by variable stars. We define our model likelihood L by 



\nL(f mt ,f h tp j ) = 



riVvar+y/jA^) 



y var i / 

1=1 



N abB 

+ E ln 

i=l 



dtFWHMidxidyi 



+ 

i=i 



d 3 N (i/jj) 



3 dtFWHMidxidyi 



(28) 



where / var is the fraction of variable stars relative to some 
fiducial model expectation number iVvar, fj and ipj are the 
lens fraction and mass function for component j, n c is the 
number of lens components and iVobs the number of observed 
events. For the disc and bulge components fj and ipj are 
both fixed, with fj = 1 and ipj given by equation (p7[), 
whilst for the Galaxy and M31 haloes ipj oc m~ 1 5(m — nij), 



as in equation (]26|), and fj and rrij axe free parameters. We 
define h with respect to the halo density normalizations in 
Table §. 
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The resolution of our simulation is insufficient to eval- 
uate reliably the third derivatives in equation (p§|), so we 
decouple the timescale and spatial distributions by comput- 
ing (dN/dtFWHM){d 2 N/dxdy) instead of d 3 N / dt^^wMdxdy 
within our fields. By averaging over spatial variations in the 
timescale distribution we are ignoring correlations which 
could provide us with further discriminatory information. 
However, in the limit of infinite data and perfect measure- 
ments we are still able to recover precisely the underlying pa- 
rameters because the average event duration is known with 
infinite precision. 

We assume that the distribution of variable stars traces 
the M31 surface brightness. In reality variable stars will be 
harder to detect in regions of higher surface brightness, so 
our idealized distribution is somewhat more concentrated 
than we should expect for a real experiment. We assume 
the timescale distribution of detectable variables is log- 
normal, with a mean and dispersion (lnfFWHivi} = 2 and 
(j(lntpwHM) = 0.5 (where £fwhm is expressed in days). 
Their timescales are therefore assumed to be typical of a 
wide range of lens masses (see Figure ^ and are thus least 
helpful as regards discrimination between lensing events and 
variable stars. 

To test the likelihood estimator we generate data- 
set realizations and compute their likelihood over a five- 
dimensional grid of models spanning a range of MACHO 
masses and variable star and MACHO fractions. For the grid 
sampling we assume uniform priors in the variable star and 
MACHO fractions and logarithmic priors for the MACHO 
masses. Since the events in the inner 5 arcmin of the M31 
disc are predominately due to stellar lenses (mostly bulge 
self-lensing) we count only events occurring outside of this 
region. 

6.2 First-season expectations 

Figure [h] shows the degree to which the MACHO param- 
eters can be recovered after one season in the optimal case 
where the data-set contains no variable stars. For the real- 
ization we have adopted a MACHO fraction of 0.25 and mass 
of 0.5 Mq for both the Galaxy and M31 haloes, and have 
set / V ar = 0. The MACHO parameters correspond to those 
preferred by the most recent analyses of the EROS an d MA- 
CHO teams flLasserre et al. 1999 ; [Alcock et al. 2000| ). Each 
panel in Figure [To| represents a two-dimensional projection 
of the five-dimensional likelihood, in which each point on the 
two-dimensional plane is a summation of likelihoods over the 
remaining three dimensions. Contours are constructed about 
the two-dimensional maximum likelihood solution which en- 
close a given fraction of the total likelihood over the plane. 
The contours shown enclose 34% (solid line), 68% (dashed 
line), 90% (dot-dashed line), 95% (dotted line) and 99% 
(triple dot-dashed line) of the total likelihood. The star in 
each plane shows the input values for the realization. 

The four panels in Figure |l^ depict the likelihood planes 
for M31 MACHO fraction and mass (top left), Galaxy MA- 
CHO fraction and mass (top right), M31 and Galaxy MA- 
CHO fractions (bottom left) and M31 MACHO and variable 
star fractions (bottom right). From the top- left panel we see 
that, after just one season, useful constraints are already 
possible for M31 parameters. In this realization the 90% con- 
fidence level spans around two orders of magnitude in MA- 



CHO mass (~ 0.05 — 10 Mq) and an order of magnitude in 
halo fraction (~ 0.1 — 1.1). The brown-dwarf regime is mostly 
excluded. In the upper-right panel we see that the Galaxy 
MACHO parameters are ill-defined after one season. This is 
unsurprising since Galaxy MACHOs are out-numbered two 
to one by M31 MACHOs and they have no signature compa- 
rable to the near-far asymmetry of their M31 counterparts. 
The panel shows a suggestive spike in the likelihood con- 
tours occurring at about the right mass range, though the 
contours marginally prefer a Galaxy halo with no MACHO 
component. The one firm conclusion that can be drawn is 
that a substantial contribution of low- mass lenses is strongly 
disfavoured by the data. The strongest constraints occur at 
~ 0.003 Mq, where the expected number of events peaks 
for a given fractional contribution. The likelihood estima- 
tor indicates that 0.003 Mq lenses contribute no more than 
~ 5% of the Galactic dark matter with 90% confidence. In 
the lower-left and lower-right panels of Figure |Io| we see the 
trade-off between M31 and Galaxy MACHO fractions and 
between M31 MACHO and variable star fractions, respec- 
tively. The lower-left panel indicates that a scenario in which 
there are no MACHOs is excluded with very high confidence, 
despite the large uncertainty in the halo fraction determi- 
nations. In the lower-right panel we see that the likelihood 
estimator has correctly determined that there is little, if any, 
contamination due to variable stars, with a 90% confidence 
upper limit of / va r < 0.03. 

In Figure 111] we show the results for a simulation over 
one season in which there are no microlensing events, only 
variable stars. We adopt A var = 100 and / var = 1 within 
the INT WFC fields. It is important to establish whether, 
in the event of there being no MACHOs, our likelihood es- 
timator is able to correctly determine a null result even if 
a significant number of variable stars pass the microlensing 
selection criteria. The four panels in Figure |ll| indicate that 
our estimator has been very successful as regards the M31 
MACHO contribution. The M31 MACHO fraction is con- 
strained with 90% confidence to be below 0.2 for lenses in 
the mass range 0.001 - 0.1 Mq and below 0.4 for MACHOs 
up to a few Solar masses. This despite a rate in variable stars 
comparable to full haloes of MACHOs. In the upper-right 
panel we see that there is considerable uncertainty in the 
Galaxy MACHO parameters, though interesting upper lim- 
its on the halo fraction are obtained for lenses in the mass 
range 0.03 — 0.1 Mq. In the lower-left panel we see that a 
non-zero MACHO contribution is preferred though the con- 
tours are consistent with the input model at about the 70% 
confidence level. In the lower-right panel we see that the esti- 
mator is able to constrain the number of variables to within 
±30% of the input value. Thus our likelihood estimator has 
provided us with not just an estimate of the MACHO pa- 
rameters but also an estimate of the level of contamination 
in the data-set. This estimate is completely independent of 
(and thus does not rely upon) additional information one 
might obtain from colour changes or asymmetry in the light- 
curves of individual events, or from follow-up observations. 

6.3 Evolution of parameter estimation 

Figure ^ shows another first-season simulation in which we 
adopt the same MACHO parameters as in Figure [l(^ but this 
time we also take A var = 100 and / var = 1. The contours in 
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Figure 10. Maximum likelihood recovery of MACHO parameters for a simulated data-set after one season. The input parameters 
for both Galaxy and M31 MACHO populations are 0.5 Mq for the lens mass and 0.25 for their fractional contribution. Here 
we assume there is no contamination to the data-set from variable stars. The four panels are two-dimensional projections of the 
five-dimensional likelihood space. The contours in each plane enclose 34% (solid line), 68% (dashed line), 90% (dot-dashed line), 
95% (dotted line) and 99% (triple dot-dashed line) of the total likelihood assuming a linear prior in the MACHO and variable star 
fractions and a logarithmic prior in the MACHO masses. The stars denote the input parameters. The four panels represent the 
likelihood in the planes of M31 MACHO fraction and mass (top left); Galaxy MACHO fraction and mass (top right); M31 and 
Galaxy MACHO fractions (bottom left); and M31 MACHO and variable star fractions (bottom right). The variable star fraction is 
measured relative to a rate of 100 events per year in the two INT fields. 



the plane of M31 MACHO mass and fraction appear largely 
unaffected by the presence of significant variable star con- 
tamination, and qualitatively resemble those in Figure |Io[ 
There is no evidence of estimator bias due to the presence 
of variables, which for our realization out-number the MA- 
CHOs from both haloes combined. However the Galaxy MA- 
CHO parameter estimation is clearly led astray by the pres- 
ence of variables, with upper limits on halo fraction pos- 
sible for only a narrow range of lens masses. The estima- 
tor nonetheless strongly excludes a no-MACHO hypothesis 
(lower-left panel) and provides a good estimate of variable 
star contamination levels. 

Figure [l^ shows the constraints after three seasons as- 
suming the same parameters as for Figure except that 
we have reduced the contamination level to / var = 0.3. A 
significant decrease in contamination would be expected as 
the increase in observation baseline permits the exclusion of 
a larger number of periodic variables. The constraints for 
M31 MACHO parameters have tightened up considerably, 
with a 90% confidence uncertainty of a factor four in halo 
fraction and an order of magnitude in MACHO mass. The 
constraints on Galaxy MACHO parameters have also sharp- 
ened considerably, allowing strong upper limits on the halo 
fraction to be made over a wide mass range, though the 
data in this case is consistent with a complete absence of 



Galaxy MACHOs. However, in the lower-left panel we see 
that the joint constraint on M31 and Galaxy MACHO frac- 
tion advocates a significant overall MACHO contribution. 
The lower-right panel also shows an accurate determination 
of contamination levels. 

In Figure ^ we depict constraints for ten seasons of 
data, comparable to the lifetime of current LMC surveys, 
with the variable star contamination level reduced further 
to / var = 0.1. The M31 MACHO fraction is now essentially 
specified to within about a factor of three, whilst the MA- 
CHO mass uncertainty is within an order of magnitude. We 
now also have a positive estimation of the Galaxy MACHO 
contribution and mass. The constraints on Galaxy param- 
eters are only a little worse than those for M31 after three 
seasons. The variable star contamination level is once again 
robustly determined. 

Figures [l^ to |l4| show that the likelihood estimator is 
able to distinguish clearly between microlensing events and 
our naive model for the variable star population. They also 
show that, given a lifetime comparable to the current LMC 
surveys, a sustained campaign on the INT should determine 
M31 MACHO parameters rather precisely and should also 
provide a useful estimate of Galaxy MACHO parameters. A 
more modest campaign lasting three seasons would provide 
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Figure 11. As for Figure [Io| but this time there are no MACHOs, only variable stars. The input model has a variable star fraction 
of unity. 
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Figure 12. As for Figure |l^, with the same input parameters, except that we now adopt a variable star fraction of unity rather 
than zero. 
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Figure 13. As for Figure 12 but for three seasons of data and a variable stars fraction of 0.3. 
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Figure 14. As for Figure |l^ but for ten seasons of data, comparable to the lifetime of current LMC surveys, and a variable stars 
fraction of 0.1. 
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a robust estimate of M31 MACHO parameters and useful 
constraints on the Galaxy MACHO fraction. 

Since all the above simulations assume the same halo 
fraction and MACHO mass for both galaxies, we decided 
to test whether our likelihood estimator was sensitive to 
Galaxy MACHO parameters independently of M31 MACHO 
values. We therefore ran a simulation over three seasons in 
which 30% of the M31 halo comprises 0.5 Mq lenses and 
60% of the Galaxy halo comprises 0.03 Mq lenses. The 
Galaxy MACHOs actually out-number the M31 MACHOs in 
this model. Whilst the model is somewhat contrived, and is 



already ruled out w ith high confidence (Lasserre et al. 1999 



Alcock et al. 2000), it provides a useful test case for our es- 
timator. We find that the estimator successfully resolves the 
mass scales of the two populations within 90% confidence, 
though with a slight tendency to overestimate the Galaxy 
MACHO mass and underestimate the M31 MACHO mass. 
Whilst we find a large overlap in preferred halo fraction, 
this is consistent with the sensitivity typically achieved after 
three seasons when the MACHO masses in the two galaxies 
are the same. The Galaxy MACHO parameters are much 
better defined than in Figure [l3|, though for this case the 
variable star contamination level was set to zero. 

There is one aspect, however, in which our simulation 
presents an over-optimistic picture. The success of the es- 
timator in discriminating between variable stars and mi- 
crolensing events is mostly due to the fact that our adopted 
variable star distribution is significantly more concentrated 
than the microlensing distribution of either M31 or Galaxy 
MACHOs. The assumption we have made, that their ob- 
served distribution traces the M31 surface brightness, is 
reasonable only for very bright variable phenomena which 
would be detected regardless of where it occurred in M31. 
For less prominent variables there will be a bias against their 
detection in the central regions of M31, where the surface 
brightness is high and so their contribution to the super- 
pixel flux relatively small. We might well expect a realistic 
distrib ution of variable stars to resemble that of Galaxy MA- 



CHOs because the surface density of Galaxy MACHOs does 
not vary significantly over the M31 disc, so their distribution 
would also trace the M31 surface brightness if there was no 
detection bias away from regions of high surface brightness. 
However, in the absence of a conspiracy between the flux dis- 
tribution of variable stars at peak luminosity and the flux 
distribution of microlensed sources at peak magnification, 
there should be some distinction between the spatial distri- 
butions of Galaxy MACHOs and M31 variable stars, though 
this may be only mild. In any case, even if the two distribu- 
tions are indistinguishable this should not significantly affect 
the determination of M31 MACHO parameters because the 
likelihood relies heavily on evidence of near-far asymmetry 
(which is why the likelihood contours are much better de- 
fined for M31 MACHOs than for Galaxy MACHOs). This 
cannot be replicated by variable stars. Only if several hun- 
dred variable stars passed the selection criteria every season 
would the signature of asymmetry be washed out and the 
constraints on M31 MACHO parameters severely degraded. 
Such an occurrence would warrant critical re-examination of 
the selection criteria! 



7 CONCLUSIONS 

Pixel lensing is a relatively new and powerful method to al- 
low microlensing searches to be extended to targets where 
the sources are unresolved. It heralds the possibility of de- 
tecting or constraining MACHO populations in external 
galaxies. Though pixel lensing is hampered by changes in 
observing conditions, which introduce spurious variations in 
detected pixel flux, techniques have been developed which 
minimize these variations to a level where genuine microlens- 
ing signals can be detected. 

POINT- AGAPE and another team (MEGA) have em- 
barked on a major joint observing programme using the 
Isaac Newton Telescope (INT) to monitor unresolved stars 
in M31 for evidence of pixel lensing due to MACHOs ei- 
ther in the Galaxy or M31 itself. Two techniques, the Pixel 
Method and difference imaging, are available to minimize 
flux variations induced by the changing observing condi- 
tions. In this paper we have assessed the extent to which the 
Pixel Method allows us to determine the mass and fractional 
contribution of MACHOs in both M31 and the Galaxy from 
pixel-lensing observables. Our assessment takes account of 
realistic variations in observing conditions, due to changes 
in seeing and sky background, together with irregular sam- 
pling. 

Pixel lensing observables differ from those in classical 
microlensing, where one targets resolved stellar fields, in that 
one is generally unable to measure the Einstein radius cross- 
ing time, t e , of an event. The fact that the source stars are 
resolved only whilst they are lensed means that one is unable 
to determine their baseline luminosity, so neither the magni- 
fication nor the total duration of the event can be measured. 
As an alternative to t e one may measure the full-width half- 
maximum timescale, £fwhm , directly from the light-curve. 
However, this provides only a lower limit to the underly- 
ing event duration. Fortunately, M31 provides a signature 
which permits an unambiguous determination of whe ther or 
not M ACHOs reside in its halo: near-far asymmetry ( protts 



1992). If M31 is embedded in a dark spheroidal halo of MA- 



CHOs the high disc inclination should provide a measurable 
gradient in the observed pixel lensing rate. The strength of 
this signature depends both on the mass and fractional con- 
tribution of MACHOs in M31, as well as the level of "con- 
tamination" by variable stars, M31 stellar lensing events and 
foreground Galaxy MACHOs. 

We have employed detailed Monte-Carlo simulations to 
estimate the timescale and spatial distributions of MACHOs 
in both our Galaxy and M31 for spherically-symmetric near- 
isothermal halo models. We also model the lensing contribu- 
tion due to disc and bulge self-lensing. The expected number 
of M31 MACHOs for our two INT fields peaks at about 100 
events for ~ 0.01 M MACHOs, the Galaxy MACHO con- 
tribution being about half as large. For a given mass and 
halo fraction we expect to detect about an order of magni- 
tude more events than current conventional surveys target- 
ing the LMC. 

The timescale distributions for Galaxy and M31 MA- 
CHOs are practically identical because of the symmetry 
of the microlensing geometry. Our simulations also confirm 
that £fwhm is less strongly correlated with lens mass than 
t e - For our sampling we find that, empirically, ^fwhm) oc 



1/2 



i 1 ' 4 for lens mass m. Sampling introduces a sig- 
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nificant bias in the duration of detected events with respect 
to the underlying average for very massive and very light 
MACHOs. 

Our simulations clearly show the near-far asymmetry in 
the M31 MACHO spatial distribution. However, the pres- 
ence of the foreground Galaxy MACHOs makes its mea- 
surement more difficult. We also find that the distribution 
of very massive MACHOs is noticeably more centrally con- 
centrated than that of less massive lenses. Stellar self-lensing 
events are found to be mostly confined to within the inner 
5 arcmin of the M31 disc, and are mostly due to bulge self- 
lensing. Their tight spatial concentration means that they 
do not pose a serious contamination problem for analysis of 
the Galaxy and M31 MACHO populations. 

We have constructed a maximum likelihood estimator 
which uses timescale and position observables to simultane- 
ously constrain the MACHO mass and halo fraction of both 
M31 and the Galaxy. The statistic is devised to be robust to 
data-set contamination by variable stars. We find that M31 
MACHO parameters can be reliably constrained by pixel 
lensing. For simulated INT data-sets we find pixel-lensing 
constraints on the M31 halo to be comparable to those ob- 
tained for the Galaxy halo by the conventional microlensing 
surveys. Even with severe contamination from variable stars 
the M31 MACHO parameters are well determined within 
three years. In particular, if there are few MACHOs in M31 
this should become apparent after just one season of data 
collection, even if as many as a hundred variable stars pass 
the microlensing selection criteria, because of the absence of 
near-far asymmetry. Pixel lensing is less sensitive to Galaxy 
MACHO parameters. Our simulations indicate that we re- 
quire at least three times as much observing time in order 
to produce comparable constraints on Galaxy MACHO pa- 
rameters. If the spatial distribution of variable stars closely 
follows that of Galaxy MACHOs, then it may become very 
difficult to reliably constrain Galaxy MACHO parameters. 

The work presented here clearly demonstrates that a 
vigorous monitoring campaign on a 2m class telescope with 
a wide-field camera can identify and characterize MACHOs 
in M31. We now have the opportunity to unambiguously 
establish the existence or absence of MACHOs in an exter- 
nal galaxy. The advantage of targeting M31 over our own 
Galaxy is that we have many lines of sight through the halo 
of M31 and a clear signature with which to distinguish M31 
MACHOs from stellar self-lensing, the primary source of sys- 
tematic uncertainty for Galaxy halo microlensing surveys. 
M31 therefore represents one of the most promising lines of 
sight for MACHO studies. 
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